Dephasing of a particle in a dissipative environment 



Doron Cohen and Baruch Horovitz 

Department of Physics, Ben-Gurion University, Beer-Sheva 84105, Israel 

Abstract. The motion of a particle in a ring of length L is influenced by 
a dirty metal environment whose fluctuations are characterized by a short 
correlation distance £ <^ L. We analyze the induced decoherence process, and 
compare the results with those obtained in the opposing Caldeira-Leggett limit 
> L). A proper definition of the dephasing factor that does not depend on a 
vague scmiclassical picture is employed. Some recent Monte-Carlo results about 
the effect of finite temperatures on "mass renormalization" in this system are 
illuminated. 



1. Introduction 

What is the dephasing of a particle that has an energy E if it is subject to a fluctuating 
environment that has a temperature T? In particular what is the dephasing close to 
equilibrium (E ~ T), and what happens in the limit T — ► 0? This question has 
fascinated the mesoscopic community during the last two decades [TJ HJ [31 IH O O [7] . 
Our purpose is to study this question within the framework of linear response 
theory for a general characterization of the environment. In the Caldeira-Leggett 
(CL) framework [5J 10 the effect of the environment is characterized by a friction 
coefficient r\ and by a temperature T. But more generally []Jjl [Til [H] it has been 
emphasized that the proper way to characterize the environment is by its form factor 
S(q,u>). The form factor contains information on both the temporal and the spatial 
aspects of the fluctuations, and in particular one can extract from it not only T and rj, 
but also the spatial correlations. Typically (but not always) these spatial correlations 
can be characterized by a correlation distance I. 

So now we ask the refined question: Given S(q,Lj), what is the dephasing? But 
first we have to say what do we mean by dephasing. In Ref. |10l ITT) the CL approach 
has been generalized. Namely, it has been realized that an environment with a given 
S(q, u>) can be modeled as a set of Harmonic oscillators. Then it is possible to apply 
the Feynman- Vernon formalism in order to trace them out. Using a semiclassical point 
of view the propagator is expressed as a sum over pairs of classical trajectories. One 
observes that after time t the interference contribution (from the off diagonal terms in 
the double sum) is suppressed by a factor P v . This factor is interpreted as a dephasing 
factor, and its expression can be cast into the form P v = exp(—F(t)) with 



t. We emphasize here the Ohmic CL model, which is of relevance in the present context as a limiting 
case for a dirty metal environment. Obviously in general one may consider non-Ohmic models, where 
memory kernels are involved while the notion of a friction constant rj becomes ill defined. 




(1) 
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where the dq integration measure depends of the dimensionality. In the semiclassical 
treatment S(q, to) would be the symmetrized form factor of the environment and 
P(q,uj;t) would be the classical symmetric power spectrum of the motion. There 
are some subtleties in the definition of P{q, u>; t) that we are going to discuss later on. 
In particular we note that P(q, w; t) may have weak dependence on t because a finite 
time interval is considered. 

It has been further argued in Ref. 12 that due to inherent limitations of the 
semiclassical (stationary phase) approximation the physically correct procedure is to 
use the non-symmetrized quantum versions of S(q,u>) and P(q,uj;t). This point has 
been further discussed in [13]. However, a proper derivation of Eq.([T]), that does not 
rely on the semiclassical framework, has not been introduced. One objective of the 
present paper is to extend the derivation of the above formula beyond the semiclassical 
context. 

It is important to realize that Eq. |T]) is capable of reproducing all the established 
results about dephasing in normal metals, including the high temperature cx T 
dependence of the dephasing time, the low temperature cx T 3 / 2 dependence in the 
case of a diffusive particle, and the cx T 2 dependence in the ballistic regime. At 
finite temperatures there is a finite time scale h/T that allows the the approximation 
F(t) w T v t and hence the notion of dephasing rate Y v is well defined. In the limit 
T — > we always have T v — > 0. This, however, does not exclude sub-exponential 
(power law) decay of P v . 

Indeed it is well known [6] that for a Brownian particle with CL environment the 
function F{t) grows as log(t) at zero temperature implying sub-exponential dephasing 
at T = 0. It is a common miss-conception that Eq.JT]) with non-symmetrized spectral 
functions fails to reproduce this log(t) spreading. We shall dwell on this point later 
on in this paper. Furthermore we shall study whether similar sub-diffusive behavior 
can be found for general S(q,uj). 

During the last decade the study of a particle in a ring, has become a paradigm 
for the study of ground state anomalies. [HI [HI EE HZl HB] . Besides being a prototype 
problem that possibly can be realized as a mesoscopic electronic device, it is also of 
relevance to experiments with particles that are trapped above an "atom chip" device 
[19l [20[ [21] , where noise is induced by nearby metal surfaces. A significant progress 
has been achieved in analyzing the equilibrium properties of this prototype system, 
in particular the dependence of the ground state energy on the Aharonov Bohm flux 
through the ring. The derivations of the dephasing factor using Eq.([T]) for the ring 
problem is a major objective of the present paper. In this context there is a growing 
understanding that the study of dephasing is intimately connected with the study of 
mass-renormalization at low temperatures. We believe that our results shed new light 
on some recent findings [22] that have been obtained using Monte-Carlo data for the 
temperature dependence of the mass-renormalization effect. 

Outline: In sections 2-3 we characterize the environment by the power spectrum 
of its fluctuations, and then in Sections 4-6 we derive the formula for the dephasing 
factor. This formula is applied in Sections 7-8 to the calculation of the dephasing 
of a particle of mass M in a ring of length L. The results depends crucially on the 
correlation distance I of the fluctuating environment. They shed light on some new 
findings regarding mass-renormalization in this system as explained in Section 9. It 
is conjectured that the mass-renormalization effect involves a measure for coherence. 
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Some further discussion of the theoretical framework is presented in Section 10. 



2. The characterization of a fluctuating field 

We regard the environment as a fluctuating field U(x,t). See appendix A for 
more details on its Hamiltonian modeling. The fluctuations of the environment are 
characterize by a form factor: 

S(q,u) = J J [(W(jc,t)«(0,0))] % lu}t ~ lqx dtdx (2) 

where the expectation value assumes that the bath is in a stationary state of its 
unperturbed Hamiltonian. The force operator is formally defined as T = —U'{x,t), 
where the derivative is taken with respect to x. The force-force correlation function is 
obtained via double differentiation of the correlation function. In particular the local 
power spectrum of the fluctuating force is 

%;) = J ^<?S{q,w) (3) 
An Ohmic environment is characterized by 

SowcM = 1 _ 2 e Xr fOT M < -c (4) 

The friction coefficient characterizes the response of the environment to an x variation 
("force proportional to velocity"). Accordingly it is given by the Kubo formula: 



1 



2ftid 



S(uj=Q) 



V = Km — S(co)-S(-lo) = -L_ ± (5) 



2T 



For a strictly Ohmic bath the friction coefficient is frequency independent and the first 
equality holds for any ui < ui c (no need to take a limit). The generalization of the above 
to 3 dimensions is straight forward. The position coordinate becomes x = (x,y,z) 
and accordingly q should be replaced by q — (q x ,q y ,q z ), with integration measure 
d 3 q/(2ir) 3 . In the definition of S(oj) the q 2 should be replaced by q 2 . or optionally by 
<7 2 /3. The simplest type of environment is known as the CL model, where the particle 
interacts with long wavelength modes. The associated form factor is 

^CaldeiraLeggett 

(g,w) = S ohmic M x 3 V > 2 KH > (6) 

Another case of interest is the interaction with a dirty metal (Appendix B) for which 

4tt£ 3 1 
5 DirtyM etai(9,w) w S ohmic (w) x — — for \q\ < - (7) 

q z I 

where the friction coefficient can be expressed in terms of the conductivity: 
e 2 1 

" = 7 X 4^ (8) 
In the latter context it is customary to define a dimensionless friction coefficient as 
follows: 

1 2 ^ ^ , , 

a = 2^ = 8^1 = 8(k F £) 2 {) 
The motion of a classical Brownian particle of mass M under the influence of such 
fluctuating environment is characterized by a damping rate 

V 27TQ! 

7 = = MP (10) 
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3. The fluctuations within a ring 

In the present paper we consider a particle in a ring of radius L/(2tt). We assume that 
S(q,uj) — S ohln i C (u))w(q) is factorisable, as in the examples of Eqs. (|6|7p . Since the 
motion is confined to one dimension it is natural to expand the spatial correlations 
along the rings in Fourier series: 

y^tf( g )e*-[«C*>-«<«0] = jr w m e im ^- e ^ (11) 



m— — oo 



Accordingly, using Eq.© with x = (L/2tv)9 we get 

OO 

S(q,Lu) = 5 ohmic (w) x 2J w m 2ir8(q- q m ) (12) 

m— — oo 

where the discrete modes are 

2ir 

q rn = -y-m, m = 0,±1,±2, ... (13) 

By convention we want r\ to be the friction coefficient. Therefore S(lu) as defined by 
Eq.Q should equal S ohlnic (ui) of Eq.Q. This implies the following sum rule: 

oo 

W m q 2 m = l (14) 

m— — oo 

In general we have ~ (L/£) fluctuating modes, each has the weight w n ~ £ 3 /L. In 
Appendix C we show that for a CL bath we have only one fluctuating mode (|m| = 1) 
with 

w ™=\{i) 2 (15) 

while in the case of a Dirty metal with short range correlated fluctuations we have 
A4 = (L/(2ir))/£ 3> 1 fluctuating modes with weights 

I 2 1 { M\ 
w m fa — x — In — - for \m\ < M (16) 

2ir M \\ni\J 

In both cases we ignore the m=0 mode for a reason which is explained in the next 
section. It is important to realize that the CL model can formally be regarded as a 
special limit of a dirty metal environment with I ^> L. In the latter case the weight 
of the | to | > 1 modes is smaller by powers of L/£ (Appendix C). 



4. The dephasing factor 



The dephasing factor P v is a number within [0, 1] that characterizes the suppression 
of coherence. Its popular definition is based on a semiclassical picture. Using the 
Feynman- Vernon formalism the propagator is expressed as a sum over pairs of classical 
trajectories. One observes that after time t the interference contribution (from the off 
diagonal terms in the double sum) is suppressed by a factor 



U[x A ] X 



U[x B ] X 



-S N [x J 



(17) 
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where x is the preparation of the batrjj| In order not to complicate the notations, the 
canonical average over \ states is implicit. The unitary operator U[x] generates the 
evolution of the bath given that the particle goes along the trajectory x(t). The action 
Sn[x a , x b ] is a double time integral. Using manipulation as in Ref. |ll[[T2"] one obtains 
Eq.JT]) with the symmetrized version of S(q,ui), and the symmetric classical version 
of P(q,ui). This semiclassical expression is definitely wrong for short range scattering 
at low temperatures |12) . because it does not reflect that closed channels cannot be 
excited. This problem with the semiclassical (stationary phase) approximation is well 
known in the theory of inelastic scattering. One way to overcome this limitation is to 
refine the definition of the dephasing factor using a semiclassically inspired "scattering" 
point of view as in Appendix D. However it is clear that such a refinement makes the 
concept of dephasing too vague. 

We propose here a natural definition for the dephasing factor that is related to 
the purity trace(p 2 ) of the reduced probability matrix. The notion of purity is very 
old, but in recent years it has become very popular due to the interest in quantum 
computation [23] . Given that the state of the system including the environment is 
Vl/pn, where p and n label the basis states of the particle and the bath respectively, 
the purity is given by 



P v (t) = y / trace(p 2 ya ) = v /trace( ) o 2 nv ) 



J2 ^P>n>%»n^ P "n»% ln » 
p' p"n'n" 

Assuming a factorized initial preparation as in the conventional Feynman- Vernon 
formalism, we propose the loss of purity (P v < 1) as a measure for decoherence. A 
standard reservation applies: initial transients during which the system gets "dressed" 
by the environment should be ignored as these reflect renormalizations due to the 
interactions with the high frequency modes. Other choices of initial state might 
involve different transients, while the later slow approach to equilibrium should be 
independent of these transients. In any case the reasoning here is not much different 
from the usual ideology of the Fermi golden rule, which is used with similar restrictions 
to calculate transition rates between levels. 

Writing the initial preparation as ^ p °2 — S PtPo S n>no , and using leading order 
perturbation theory, we can relate to the probabilities P t (p,n\po,no) — \^ pn \ to 
have a transition from the state \po, no) to the state \p, n) after time t. The derivation 
is detailed in Appendix E. One obtains the result 

Pipit) = Pt(po,no\po,n ) + P t (pjtpo,no\po,n ) + P t (po, n^n \po, n ) (19) 

in agreement with the semiclassically inspired point of view of Appendix D. The 
notation p ^ po or n ^ no implies a summation J2 p ^ po or J2 n ^n > respectively. 
In the next section we shall discuss the actual calculation of Pt(p, n\po, no), using 
the Fermi-golden-rule (FGR). Thus we deduce that within the FGR framework, 
the purity is simply the probability that either the system or the bath do not 
make a transition. The first term in Eq. (|I9p is just the survival probability of the 
preparation P aurviva i = Pt{po, no\po, n o)- The importance of the two other terms can 



1/2 

(18) 



§ The implicit assumption of initial factorized state is of course problematic |241 1251 . In most cases 
it implies an unpleasant transient that should be ignored. We further discuss the significance of the 
long time decay later in this section after Eq. (1 1 8 I t and in the Summary 
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be demonstrated using simple examples: For an environment that consists of static 
scatterers we have -P Burvlval < 1 but P v — 1 thanks to the second term. For a particle in 
a ring that interacts with a q=0 environmental mode .P survival < 1 but P v — 1 thanks 
to the third term. Using J2 P n P{p, n\po, ^o) = 1 we obtain the optional expression 

P v = 1-P V = Yl J2 Pt(P,n\po,n ) (20) 

p^po n^no 

In the problem that we consider in this paper we can calculate P v using a dqdoj 
integral as in Eq.([T]). In many examples the uj = transitions have zero measure 
and therefore P v is practically the same as Pt{pa, «o|poj no). Otherwise one has to 
be careful in eliminating those transitions that do not contribute to the dephasing 
process. Anticipating the application of Eq.([T]) for the calculation of the dephasing 
for a particle in a ring, the integration over q becomes a summation over q m , and the 
m = component should be excluded. 



5. Dephasing at finite temperatures 

The interaction between the particle (x) and the environment can be written as in 
Appendix A: 

V = J dxp{x)U{x) (21) 

where p(x) = 5(x—x). In the Heisenberg (interaction) picture a time index is added so 
we have U{x,t) and p(x, t) . Given a preparation of the bath and of the system we can 
define S(q,u) to characterize U(x,t) and we can also define P(q,u>) to characterize 
p(x,t). The precise definition of the latter object is further discussed below. The 
survival probability of a factorized preparation is P v {t) = l—p v (t) where: 

p v (t) =11 {V{ti)V{t 1 ))dt 2 dt 1 (22) 



dt\dt 2 JJ dx\dx2 \p(x2, t 2 )U{x2, t 2 ) p(x%, ti)U{xi, t\) 

^L^-S(q,Lo) JJ dhdt 2 JJ dx,dx 2 (p(x 2 ,t 2 )p(x 1 ,t 1 ))eM*'-* 1 '>- i "&-V 

At finite temperatures, if recurrences due to finite-size quantization effect can be 
ignored, one can obtain as an approximation p v ps T ip t 1 where is called the 
dephasing rate. In the next section we discuss circumstances where such an 
approximation is not valid: the feasibility of this approximation requires neglect of 
the end-point contributions to the double time integration. By going to the variables 
(ti +t 2 )/2 and t = t 2 — t\ one obtains the following expression for the dephasing rate: 

r v - /f^(?, W )P(- g ,- W ) (23) 

The implied definition of P{q, u>) is discussed below and further refined in the next 
section. Following standard argumentation one conjectures that the long time decay 
of Ptp{t) is exponential, as in the analysis of Wigner's decay. The similarity of Eq. ([23|) 
to the semiclassical result (as discussed below Eq.([T])) is obvious. It is important 
to realize that in the present context the non-symmetrized quantum version of the 
power spectrum has emerged. Furthermore, if we want to calculate P v , and not just 
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the survival probability of the initial state, we have to be careful about the proper 
treatment of the diagonal terms as discussed in the previous section. Accordingly we 
eliminate the diagonal term from the implied definition of the power spectrum: 

f+OO 

P(q,u) = / ( e -ipW e i?i(o)^ _ (jqxy e i"r dT ( 24 ) 



We emphasize again that in a later section we are going to treat the time limits more 
carefully, where P(q,u>) will be replaced by P(q,uj;t) as in Eq.fT]). For a ballistic 
particle with mass M and momentum p — (2MB) 1 / 2 we have: 

P{q,u) = 2ttS(u; - u(q)) (25) 

where u>(q) = [(p + q) 2 — p 2 ]/ (2M). The power spectrum is illustrated in Fig. 1. The 
expectation value in Eq. (|25[) is taken for a particle with momentum p. For the ground 
state p—0 and hence u>(q) — q 2 /2M. In particular for a particle on a ring Eq. (fT2|) 
implies u){q m ) = q 2 a /2M. The ballistic case should be contrasted (see Fig. 1) with the 
power spectrum of a diffusive particle: 

*«■"> = zjIwf < 26 > 

where D is the diffusion coefficient. In the ballistic case the power spectrum is 
concentrated along oj = oj(q), while in the diffusive case it spreads over the range 
\u>\ < Dq 2 . In any case 

f duj 

I P(q,uj)—— = 1 by definition, for any q (27) 

Assuming close-to-equilibrium conditions, the expectation value in Eq. (|24| should 
reflect a thermal state with energy E ~ T. In practice one may set E ~ 0, though in 
general better to be careful about it: looking at Fig. I one can deduce that to take 
E r*> 0, in the problem that we aregoing to consider, results in an underestimation of 
the dephasing rate by a \/2 factoi[jJ. 



6. Dephasing at "zero" temperature 



The expression for T v gives manifestly zero dephasing rate in the limit of zero 
temperature, because in this limit S(q,Lj) and P(—q, —lu) have no overlap. However, 
this does not mean that P v does not decay. It still might have a sub-exponential 
decay. In order to understand this point we first discuss a simple artificial calculation 
of the double time integral 

P*>(*)= / / C(ti-t 2 )dtidt 2 (28) 
Jo Jo 

where C(r) is the symmetrized force- force correlation: its Fourier transform C(u) is 
the symmetrized version of S(lo) of Eq.(j4|). Later in this section we come back to the 
actual calculation and discuss how Eq. (|2"2")) can be treated. 

|| Assuming that E ~ we are going to explain in Section 7 that the non- negligible contribution to 
the integral comes form the range \q\ < qj< of effective modes: only the fluctuating modes in the 
rectangular region of Fig. f resonate with the particle and hence contribute. The power spectrum 
of the particle for E ~ T is shifted "upwards" in u and consequently the effective q range becomes 
larger by factor y/2 compared with the E ~ case. 
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The approximation p v (t) w T v t for the integral in Eq. (|2"5)) is based on the 
assumption that C(r) has short range correlations with non- vanishing integral. At 
finite T this assumption is indeed satisfied because S(u=0) = 2r\T is finite. But at zero 
temperature the integral over C(r) is zero. In fact at zero temperature Eq.Q implies 
that C(ui) = 2t]\l)\ and hence C(r) has power law tails — (ry/7r)/r 2 . It is important to 
realize that the T=0 behavior prevails also at finite temperatures provided T < lo c , 
and the time of interest should be smaller compared with 1/T. Under such "T=0" 
conditions C (r) can be approximated by its T=0 version. In order to see what comes 
out from Eq. (j2"8)) we observe that 



(29) 



r v (t) = - Pv (t) = I C(r)dr 



For an Ohmic bath at "zero temperature" the integral over the power law tails of C(r) 
gives r(t) oc 1/t, hence the spreading is logarithmic: 



P v {t) = —r\ ln(u) c t) + const 

7T 

It is instructive to make the same calculation in ui space. One realizes that 



p<p(t) 



sin(o;t/2) 



u/2 



(30) 



(31) 



which for C[ui) oc \w\ gives correctly the logarithmic spreading. 

Without any approximation we can generalize the above treatment so as to handle 
Eq. (|22[) . taking also into account the non-symmetrized nature of the spectral functions. 
Performing the dx\dx2 integration we obtain p v {t) — F(t) as in Eq.([l]) where 



P(q,w;t) = 



-iqx(t 2 ) e iqx(t 1 )\ 



(l) dhdt 2 



(32) 



In complete analogy with the way in which the environmental fluctuations has been 
treated, we express the correlator as a Fourier integral over P(q, u>), and then we are 
able to explicitly perform the dtidt2 integration. The outcome of this procedure allows 
to express the result as a convolution: 

2 

(*)P(g,w) (33) 



1 



sin(cjt/2) 



w/2 



An optional compact way of writing the final result is 



p v {t) 



dq 



dcj dw' ~ ~ 
T5-^-S(<?^) p (-<7>- w ) 

Z7T Z7T 



su\{{uj-uj')t/2) 



(34) 



(w-u}>)/2 

We note that in practical calculations or for aesthetic reasons it is possible to make 
the replacement 



sin(wt/2) 
lj/2 



(2/t) 



x t 



[optional 



(35) 



.(1A) 2 +< 

The more convenient Lorentzian kernel has the same normalization, the same 
width 1/t, and the same 2/uo 2 tails. It can be regarded as arising from using "soft" 
rather than "sharp" cutoff for the time integration. 
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7. Dephasing in the presence of a dirty metal "T>0" 

We turn to consider a particle of mass M in a ring of length L. We assume close-to- 
equilibrium conditions so we take the energy of the particle above the ground state as 
E ~ T. We consider in this section temperatures T that are much larger compared 
with the level spacing A ~ (ML 2 ) . This does not mean that the system is not 
coherenfl We further assume that the time of interest is much longer compared with 
the relevant dynamical time scales, and in particular compared with 1/T. With this 
assumption it is legitimate to use Eq. (|23[) to calculate the dephasing rate T v , and 
to treat the du> integration as if the levels of the ring form a continuum. With the 
substitution of Eqs. ([12l25p this leads to the following result: 

T v = w m S ahmio (-w(q m )) (36) 

m 

Of course one has to verify at the end of the calculation the self-consistency condition 
-C min{T, lo c }. This condition would be satisfied if the system-environment 
coupling is not too strong. 

A graphic illustration of the (q,u>) integration domain is presented in Fig. 1. The 
T dependence of S oSlm i C in Eq. ()36|) limits the sum to u)(q m ) < T. Taking into account 
the weight factors the effective number of fluctuating modes is 

M etf = — w min{ M, q c L, q T L } (37) 

where q c = (2Mcj c )^' 2 and qx — (2MT) 1 ' 2 . Note that if we had Fermi occupation 
it would be qr = T/vf, while for diffusive motion it would be qr = (T/D) 1 ^ 2 . The 
dephasing rate is obtained by summing over all the contributing modes. Each effective 
mode contributes 2r(T x w m to the sum. Accordingly 

T v = 2nT x w ™ ~ 2r l T x wM ett (38) 

0<|g m |<<? c ff 

where the average weight is w ~ £ 2 /A4 ~ £ 3 /L for a short range correlated dirty metal 
environment (< < L), while w ~ L 2 in the opposite CL limit (i 3> L), as implied by 
Eq. I]16p and Eq. (fT5|) respectively. If all the modes are effective we get T v = 2nT x £ 2 , 
while in the case of a CL environment we get the well known result 



T v = 2nT x L 



Caldeira-Leggett 



(39) 



For a fluctuating environment with correlation distance £, Eq. ([37)) implies a crossover 
temperature: 

T* = minj^, Wc | (40) 

For T < T* the dephasing rate depends on qx and therefore develops non-linear 
dependence on the temperature, as illustrated in Fig. 2 and further discussed below. 
Using a field theoretical approach [T5] it is argued that the renormalized value of the 
high frequency cutoff is 



= max{ 7, A } (41) 

effective 



As discussed in section 9 the coherence measure is r^/Aeff. In the regime of main interest Y v T. 
Furthermore, the relevant energy scale is not necessarily the level spacing: in the mesoscopic context 
the relevant energy scale (e.g. the "Thouless energy" ) is typically much larger and proportional to h 
in contrast to the microscopic quantization scale which is proportional to /j dlmBnalonaht y i 
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where 7 is the classical damping rate of Eq. lfTU)) . and A is the level spacing. The 
reasoning is as follows: all the higher frequencies contribute to mass renormalization 
only, and do not affect the dephasing process. In more details: The significant 
renormalization starts only below lu c where the linear \lo\ dispersion of the dissipation 
term dominates and leads to In a; terms in perturbation theory and to the need of 
either renormalization group or an equivalent variational method [18j . 

We would like to remark that if we do not apply Eq. ([4Tj) . the results that we 
derive below are affected quantitatively but not qualitatively. Substitution of Eq. (|4Tj) 
into Eq. (|4"0"|) implies that r\£ 2 < 1 and r\£ 2 > 1 define distinct regimes of behavior. For 
a dirty metal environment r\£ 2 <C 1 is equivalent to a <C 1 i.e. Uf£ ^S> 1. For a CL 
environment r\£ 2 is formally infinite, or one may say that i is effectively determined 
by the finite size L of the system. 

We come back to the dephasing rate calculation. In the case of a fluctuating 
environment with a short correlation distance £, we see that the high temperature 
(T > T*) result is: 

f (2r]£ 2 ) T if rj£ 2 > 1 

\ (2n£ 2 fl 2 T if rj£ 2 < 1 

where in the r)£ 2 <C 1 expression we have identified the effective (renormalized) cutoff 
as uj c = 7- The strong coupling result {q£ 2 3> 1) cannot be trusted because the self- 
consistency requirement (T v <C T) is not satisfied. This is not in contradiction with 
the observation that the CL result Eq. (|3T)j) is formally a special case of the strong 
coupling result with £ h- >• L. In the latter case the self consistency relation becomes 
r\L 2 <C 1 irrespective of £. 

In the low temperature regime we have (be definition) q* = qr- Consequently 
the T dependence becomes non-linear, and we get T v m n£ 3 M X / 2 T^/ 2 . The similarity 
of the latter to the familiar result for a diffusive electron is misleading. In both 
cases qr oc T 1 / 2 but for different reasons, and with a different prefactors. For sake 
of completeness we write the precise expression which is obtained for a dirty metal 
environment using Eq. ([38|) with Eq. (fT6|) and Eq.®: 



for T>T* 



(42) 



4ir 2 a \Qt£ J 



for A < T < T* 



(43) 



Since T v ~ T 3 / 2 at sufficiently low T the condition T v <C T is valid even for strong 
coupling r\£ 2 > 1. We note also that inclusion of a q = mode in l|38p would have led 
to T v ~ wqT. Hence our precise formulation in Eq. (|45[) is essential. 

The crossover from the high temperature result to the low temperature result 
is illustrated in Fig. 2. The illustration assumes r\£ 2 -C 1 which implies that the self 
consistency requirement (T v *C T) is globally satisfied. It also should be realized that 
the low temperature regime A < T < T* exists only for a "large ring" {j]L 2 3> 1), for 
which T* = 7 ^> A. Finally, for the L < £ case, the weight w±i dominates, leading to 
M e! f = 1 and T v = 2rjL 2 T as in Eq. ([39|) . which for ui c < A is consistent with < A 
for r/L 2 < 1. 



8. Dephasing in the presence of a dirty metal "T=0" 

We would like to discuss the "zero temperature" regime. One should be very careful 
in specifying the conditions of physical interest, else the problem may become trivial 
or of no experimental relevance. In what follows we assume that the dimensionless 
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coupling between the system and the environment (r\l 2 for a dirty metal, or r/L 2 for a 
CL environment) is much smaller than unity. This means that the competing energy 
scales are the level spacing A ~ 1 /(ML 2 ) and the temperature. So the simplest 
definition of zero temperature is T< A for which the system is in the ground state 
with an exponentially small probability to be found in an excited state. In this regime 
mass renormalization effect can be calculated using second-order perturbation theory, 
as in the Polaron problem, or possibly using field theoretical methods. Furthermore 
in this regime we can treat the dephasing problem using a "two level approximation" , 
which is a very well studied model [25) . 

The notion of "zero temperature" is also applicable if T ^> A provided the time of 
interest is short (t <C 1/T 1 ). In this regime the power spectrum S(q,ui) is the same as 
for T=0 within the frequency interval T <C ui <C ui c . Consequently the Ohmic temporal 
correlations are C(r) w — (t]/it)/t 2 within the time interval (l/w c ) <t< i^/T). As 
explained in a previous section such correlations may imply a logarithmic growth 
of Ptp(t). In view of the claim that the renormalized value of lo c is the damping 
rate 7, it follows that logarithmic spreading may arise only if T <C 7, which is the low 
temperature regime. 

As discussed in a previous section, Eq. (|34p gives a non-zero result for p v (t) even 
at zero temperature. For the CL model we have only q^O fluctuating modes and we 
get the expected log(t) spreading: 

/ T \ 2 



Pip(t) = — ( ^- J ln(w c i) C aldeira-Leggett (44) 

7T V 27T J L 

where lo c is the high frequency cutoff of the temporal fluctuations. More generally, for 
a particle in a ring the result can be written as a sum over all the q Fourier components: 

Mt) = ??" mln ((iA) + -(^)) (45) 

Strictly speaking for a finite system the assumption T> A always breaks down in the 
zero temperature limit. Still it is meaningful to formulate a condition for not having 
dephasing at zero temperature irrespective of the finite size effect: 

lim p v (t = 00) -C 1 (46) 

L — >oo 

Thus the question is simply whether in the continuum limit the q summation in 
Eq. (|45p converges in its lower limit. For a fluctuating environment with a finite (short) 
correlation distance I. 



p v ~ / dqr\i In 



MbJr 



<1 2 



v e In I ) (47) 



M 

We conclude that for a dirty metal environment with f «L, coherence is maintained 
if ?7€ 2 < 1, i.e. if k F i> 1. 

It is important to notice the following: in the strict Caldeira-Leggett limit (£ = 00) 
the size of the ring L can be arbitrarily large, hence the Heisenberg time 1/A becomes 
huge, and Eq. (|45|) leads to Eq. ([44"|) . which becomes an exact result. So if we consider 
a dirty metal environment with long wavelength fluctuations (£ 3> L) , it looks as if we 
are back in the "CL regime" leading to Eq. (|44"]) . But this is not quite correct unless 
we give away the weak coupling assumption r\i 2 C 1. As long as we keep a. <C 1 
(fixed) the CL result does not apply. This is because once i — > L and r/L 2 < 1 the 
quantization of the energy spectrum becomes important and the renormalized cutoff 
frequency Eq. (|4"Tj) becomes lo c ~ A instead of u c ~ 7. Accordingly, in the latter case, 
the time during which the log spreading prevails diminishes. 
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9. Mass renormalization 

It can be shown [lOj that a particle that interacts with a fluctuating 'dirty' environment 
acquires an additional inertial (polaronic) mass. However in recent works [16j H3 HH] 
the mass renormalization concept appears in a new context. The free energy F(T, $) 
of a particle in a ring is calculated, where T is the equilibrium temperature and $ is 
the Aharonov Bohm flux the through the ring. Then the coherence is characterized 
by the "curvature", which is a measure for the sensitivity to <£>. If the interaction with 
the environment is negligible the result can be written as 

8 2 F e 2 

' -f{M*L 2 T) (48) 



with the bare mass M* = M. The dependence of the curvature on T merely reflects the 
Boltzmann distribution of the energy. In the presence of coupling to the environment 
it turns out that M* > M. At T = 0, for fixed r\l 2 < 1, Monte Carlo data show 
[22] that the ratio M* /M is independent of the radius beyond a critical L c . As the 
radius becomes smaller than L c , the ratio M* /M rapidly approaches unity. In the 
regime of "large L" the mass renormalization effect diminishes with the temperature 
and depends on the scaled variable LT, while for "small L" the ratio M* /M grows 
with the temperature, and depends on the scaled variable L T. 

The natural question is whether we can shed some light on the physics behind this 
observed temperature dependence of the mass renormalization factor. In particular 
we would like to explain why in in one regime M* /M is a function of LT, while in 
another regime it is a function of L A T . Making the conjecture that the temperature 
dependence of M* /M is determined by dephasing it is natural to suggest the following 
measure of coherence: 

X(T ,L) = . Ll (49) 

Namely, it is the dephasing factor at the time t — 1/A cff , where A eff is the "relevant" 
energy scale. Equivalently the condition x <C 1 means that the energy levels near 
A cff remain sharp. The inequality x < 1/2 can serve as a practical definition for 
having coherence. It can be translated either as a condition on the temperature, 
or optionally it can be used in order to define a coherence length that depends on 
the temperature. The conjecture is that y = M*/M is a function of x. Let us 
calculate x using Eq. ([55|) . We assume i]£ 2 <C 1 but r/L 2 can be either larger or smaller 
compared to unity. This is equivalent to saying that the damping rate 7 can be either 
larger or smaller compared with the level spacing A. Using Eq. ([4T))) with Eq. ([4"Tj) this 
further implies that M ct{ is either larger or of order unity respectively. The transitions 
that are associated with the "relevant" energy levels are characterized by q ~ q oi{ and 
accordingly A cff ~ A4 ett x (ML 2 ) -1 . Using Eq. ([38]) we deduce that the result does 
not depend on M Bii but only on w, leading to 

(T t\ -MTirn (vM£ 3 LT forL»^ 

x(T,L) = rjwML T = j ^ £ <<; £ (50) 

We recall that our T v is valid at least for weak coupling rji 2 1. The scaling of T v 
with L is consistent with Monte Carlo exponents for the coherence length L ~ T _Al 
with either \i = 1 or fi = 1/4. The Monte Carlo data has not determined so far 
whether the transition between the two regimes is at L sa t or whether it is coupling 
dependent. 
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10. Summary and Discussion 

In this paper we have defined and calculated the dephasing factor P v (t) for a particle 
in a ring due to the fluctuations of a dirty metal environment. At finite temperature 
we have calculated the dephasing rate T v . Our interest was mainly in the mesoscopic 
regime A -C -C 7, where interference is important (because Y v <C 7). Unlike the 
microscopic regime (T v <C A), which is customary in atomic physics studies, here the 
question of dephasing at low temperature is tricky both conceptually and technically. 

The decoherence is induced because the system gets entangled with the 
environmental modes. It should be clear that in generic circumstances the coupling 
always induces "transitions" that lead to system-bath entanglement. Accordingly we 
have P v (t) < 1 even if "T = 0". This by itself does not mean "having dephasing": 
entanglement is also associated with the adiabatic renormalization due to the 
interaction with the high frequency modes. In order to "have dephasing" the loss 
of purity should not be just a transient: rather it should be a progressive process. 

Still even with this careful point of view, the reader may doubt whether the 
notion of "dephasing factor" is really helpful in studying dephasing. After all what 
do we "really" mean by dephasing. Possibly the "correct" procedure is to study an 
equilibrium correlation function C(t), and to ask whether it is damped in the t — > 00 
limit. In the absence of coupling to the environment the Fourier transform C(io) is 
a sum over delta functions S(u> — &l r ). Due to the coupling the deltas are broadened 
into resonances of with T r . This is true at any temperature, also at "T = 0". The 
controversy about dephasing at "T = 0" is related to the limit L — > 00. Do the 
resonances overlap in this limit? Do singular features of the uncoupled system survive? 
For sub-Ohmic bath [25] the ratio T v /A, where A is the mean level spacing, diverges 
as L — > 00. But the Ohmic case is "marginal" and within the framework of the Fermi- 
golden-rule it remains a constant a. So if this a is smaller compared with unity, we 
naively expect no dephasing at "T = 0" . 

The naive expectation of having no dephasing at "T = 0" is not without loopholes. 
One obvious loophole is the mass renormalization issue. If hypothctically the 
renormalized mass and hence the density-of-states diverge as T —> 0, it might imply 
dephasing at zero temperature. The recent studies of equilibrium properties of the 
ring problem are aimed in studying this question carefully, in a controlled way. For 
a particle that interacts with a dirty metal environment we believe, on the basis of 
P~8j [22] that the renormalized mass at zero temperature is finite. So in the case of a 
dirty metal environment there is no indication for "dephasing at T=0" . 

Still one would like to know what happens at low but finite temperature. As we 
said previously, no doubt that study of equilibrium properties is conceptually the best 
procedure. Still, we also want to physically understand the results. Here we come 
back to the 'dephasing factor' notion. In spite of the problems which are associated 
with this concept, we believe that it is powerful enough to shed light on the physics 
of dephasing. Our aim in this paper was to maximally exploit this notion, within the 
Fermi golden rule picture, in order to demonstrate that it captures the correct physics 
of all the established results regarding dephasing. In particular it has provided an an 
explanation for the T dependence of the mass renormalization effect, and under what 
conditions the spatial aspect of the fluctuations is capable of suppressing the "T=0" 
power law decay of coherence. 
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Appendix A. Modeling of a fluctuating environment 

It is customary to write the system-environment Hamiltonian as 

H t o t ai = H Bye (x,p) + K nt (x,Q a )+H exxv (Q a ,P a ) (A.l) 

where (x,p) are the canonical coordinates of the particle, and (Q ai P a ) are the 
environmental degrees of freedom. In the case of an interaction of a particle with 
a dirty metal environment (in 3 dimensions): 

V m = J d 3 xp(x) J d 3 x' ^ X J^ = J d 3 x'p(x)U{x) (A.2) 

where the electronic density n(x) can be expressed as a function of their coordinates, 
while p(x) = S(x — x) is a particle related field operator. 

In order to allow a Feynman- Vernon treatment it is more convenient to regard 
U(x) as arising from the interaction with a bath of harmonic oscillators [ID] . Each 
harmonic oscillator is a scatterer which is characterized by its location x a and 
its natural frequency Lu a . The interaction of the particle with the a scatterer is 
Q a u(x — x a ), so we write (in one dimension): 

K„t = c a Q a u(x - x a ) = dxp(x)U(x) (A.3) 

where the c a are coupling constants, and p(x) = 5{x — x). In the Heisenberg 
(interaction) picture a time index is added so we have p(x, t) and U(x, t). Accordingly 
the fluctuating filed is 



U(x,t) = ^2c a Q a (t)u(x - x a ) (A.4) 

ct 

As explained in Ref.flO] it is possible to postulate the interaction u(r), and the 
distribution of the parameters (x a , w a , c a ), such as to obtain a fluctuating field 
with a physically desired S(q,u>). This type of modeling is equivalent to the field- 
theoretical assumption of having Gaussian fluctuations, and accordingly a linear 
response treatment of the environment becomes exact. 



Appendix B. The fluctuations of a dirty metal 



For a metal we can use FD realation in order to relate the spati-temporal power 
spectrum S(q, u>) to the conductivity: 



S(q,u,) 



Aire 



-Im 



-1 



e(q,uj) 



21 1 



I _ e -hu/T 



where 



and 



e(q,u) = 1 



Ana 



-iuj + Dq 2 



Im 



_e{q,u) 



Airaui 



(Dq 2 + Attij) 2 + uj 2 



Thus we get 



e 2 1 2hw 
7^1-e-^/ T 



for \u\ < lo c , \q\ < 



1 



(B.l) 
(B.2) 
(B.3) 
(B.4) 
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The ohmic behavior is cut-off by the Drude collision frequency io c , and the elastic 
mean free path is £ = vp/oj c , where Vp is the Fermi velocity. The expression for 77 in 
Eq.© is obtained from the Kubo formula Eq.([5]). 

Appendix C. Calculation of Fourier components 

We are interested only in fluctuations within the ring. Therefore we have to calculate 
the Fourier components of the correlator (M{x(6 2 ),t2)U{x{6\),t\)). For the CL model 
we use the integral 

d 3 q e lqR = const - = const + (Jl\ cos(0a-0i) (C.l) 

where R = R(0 2 )-R(0i) so that \R\ = |2sin((0 2 -0i)/2)|[Z,/27r]. For a dirty metal 
with fluctuations within q < Xjl we have 



d 3 Q 47r£ 3 ctq , R _ £ 3 = p2 



(2tt) 3 q 2 ^R2 + p 



OO 

a o+^2 a m cos(m(0 2 -0i)) 

m— 1 



(C.2) 



Our a m=j LQ are half the "convention" in Ref . |16j . From the Fourier transform relation 

it follows that J2m=o a ™ = ^ anc ^ we a ^ so nave the sum rule J2m=i a m'm 2 = M 2 
where M. = (L/(2tv))/£. Disregarding the m = Fourier component the following 
approximation can be obtained [TBI IT?] for M, ^> 1 

' ln(— ) forO<m<M (C.3) 



nAi \ m J 

From Eq. (fTT|l it follows that w m — £ 2 a m /2. We conclude that the particle in the ring 
experiences white fluctuating field that is characterized by a correlation distance £. 
The fluctuating field can be re-interpreted as arising from a short range interaction 
u(r) with uniformly distributed set of scatterers as in Eq. (|A.4| ). 

In the other extreme case of CL-like environment {£ 3> V) the fluctuations of 
the higher (to > 1) modes are negligible compared with the fluctuations of the lowest 
(7n=l) mode. Accordingly we say that the number of effective modes is M = 1. Using 
the ad- hock notation A4 = (L/(2ir))/£ <C 1, the sum rule which is based on Eq. (IC.2|) 
implies that a\ — M 2 while a m> i have higher powers of M. 

Appendix D. The dephasing factor - semiclassical perspective 

For short range scattering, if x A and x B of Eq. ([T7|) are well separated, and hence 
interact with different sets of oscillators, we can argue within the semiclassical 
framework that P v is the probability to induce an excitation in the bath (i.e. "to 
leave a trace in the environment"). The argument is elaborated in Appendix C of 
|llj . This argument fails if the interfering states are not well separated in space, 
but rather interact with the same oscillators. For this reason the definition of the 
dephasing factor has to be refined. One possibility is to adopt a "scattering" point 
of view, hence treating correctly closed channels. Using sloppy notations the idea is 
to define the dephasing factor in analogy with Eq. fTTj) as P v = \U[ij) A ]x U[if) B ]x 



where ip A and i(j b are ingoing states of the system. This way of writing is suggestive 
rather than exact. Referring to a superposition preparation of the ring, where ip 
and tp B are momentum eigenstates, it is clear that P v is not necessarily the same as 
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the probability to induce an excitation in the bath. This is because U[\jj A ] and U[4> B ] 
involve the excitation of the same oscillators, rather than different sets of oscillators. If 
the factorized preparation is \pqu ), then we write the evolved state in the interaction 
picture after time t as |(po^o)t)- If we have initially a superposition + \p2), the 
evolved state would be 

I*) = bi)®lx (1) > + |P2>«)|X (2) ) + InelasticPart (D.l) 

where the so-called relative states of the bath are 

x « = ( Pl n\( Pl n ) t ) (D.2) 

X ( n ] = (P2n\( P2 no)t) (D.3) 

The dephasing factor is 

P v = \(x (1) \x (2) )\=Pt(p ,n \p ,n )+ xP'xP (D.4) 

n(^no) 

where we assume pi ~ p 2 ~ Po- It is not difficult to realize that the same 
approximation implies that the second term equals Ptipoi n^no\poi n o) m agreement 
with Eq. (|19[) . With some further argumentation we can justify the the second term 
in Eq.(fT^|) as well. We note that this derivation parallels the semiclassical treatment 
in Appendix D of [11], where P(q,oj) is defined as the difference P\\(q,uj) — P±(q,oj). 



Appendix E. The purity based definition of the dephasing factor 



In this appendix we explain the derivation of Eq.(|19p from Eq. (fT5]) . The zero 
order term in Eq. (fT5]) is the p'=p"=po,n'=n"=no term. It is equal to P5 2 where 
Po = Ptipdi n o|POj n o)- There are four sets of first order terms: The sum of the 
p'=p"=p ,n / =n Q ,n"^n terms is P x p sys where p sys = Pt(po,n^n \p a ,n Q ). Here 
n =/= no implies a summation ^ n _i no - The sum of the p' =p" =po, n'^no, n"=no terms 
is the same. There are two other sets, with either p 1 ^ p or p" ^ p , that give 
each P x p onv , where p Bny = Pt(p^Po, no\po, no). Summing over all the leading order 
contributions we get 



P v = 



Po 2 + 2P Q x p sys + 2P x p onv + 0( P 2 ) 



1/2 



(E.l) 



leading to P v m P + p eyB + p cny which is Eq. JT 
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Fig.l: The (q,u)) plane. The power spectrum S(q,u>) is distributed prc-dominantly within the 
rectangular area q < (l/£), that has a high frequency absorption cutoff ui c , and a lower emission 
cutoff. The emission cutoff T in this illustration reflects an assumption of having T < ui c , otherwise 
it would be equal ui c too. The power spectrum P(q, ui) which is associated with the ballistic motion 
(left panel) or with the diffusive motion (right panel) of the particle is illustrated by the dark region. 
In both cases the energy E of the particle implies a frequency cutoff, which is analogous to T. Close 
to equilibrium one should take E ~ T, but without much error we take for low temperatures ballistic 
motion E ~ 0, which is also illustrated in the left panel. 



A 




Fig. 2: Illustration of the dependence of the dephasing rate T on the temperature T. The 
dephasing rate is well defined for t > (f/T), and hence the self consistency requirement is T < T. 
This condition is demonstrated by a comparison with the dashed line. The illustration reflects an 
assumption of having r)£ 2 <C f , and therefore the crossover temperature T* is equal to the damping 
rate 7 = rj/M. The illustration further reflects an assumption of "large ring" (r)L 2 S> f ) for which 
7 > A, else the low temperature regime (A «T« T*) disappears. For extremely low temperatures, 
such that T is smaller compared with the spacing A = 1/(ML 2 ), the probability to excite the system 
is exponentially small and the familiar two-level modeling becomes applicable. 



